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Abstract 

We consider a quantization of relativistic wave equations which al- 
lows to treat quantum fields together with interacting particles at a finite 
time. We discuss also a dissipative interaction with the environment. We 
introduce a stochastic wave function whose dynamics is determined by 
a non-linear Schrodinger-type evolution equation in an additional time 
parameter. The correct classical limit requires the proper time interpre- 
tation of the time parameter. An average over the proper time leads to 
the conventional quantum field theory of particles which are free at an in- 
finite space separation. We consider models with scalar and vector fields 
on a pseudoriemannian manifold. A quantization of the Einstein gravity 
in this approach is briefly discussed. 

PACS NUMBERS: 03.70.+k,02.50.Ga,04.60.-m 

I Introduction 

The classical non-relativistic mechanics describes system's history x t as a func- 
tion of time ( a succession of events). In a classical relativistic theory time 
and space should be treated on an equal footing. If so then we should con- 
sider world lines (histories) x /x ( > ) rather than events. In quantum mechanics 
the wave function tp t (x) gives an amplitude of probability of detecting a parti- 
cle at the position x. Only in the classical limit |i/>t(x)| 2 w |7/>(x t )| 2 we regain 
the time-dependent trajectory. By an analogy to the non-relativistic quantum 
mechanics in an explicitly relativistic treatment we would require that in the 
classical limit the probability of an occurence of the history x^(-) in a state ip 
should be \ip(X)j,(-))\ . This probability amplitude can vary depending on a cer- 
tain universal time r which should be a Lorentz scalar. In classical relativistic 
physics there is a good candidate for this universal time (the proper time) which 
describes e.g. the periods of oscillations of a monochromatic light emitted by 
an atom or the frequencies of crystal vibrations in their rest frames. 

The relativistic quantum mechanics is usually rejected in favor of the quan- 
tum field theory (QFT). There is at least one good reason for that: the need 
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to describe the processes of particle creation and annihilation. The QFT suc- 
cessfully and consistently treats the scattering processes with particle creation 
and annihilation. However, in order to describe bound states (e.g. the Lamb 
shift in the hydrogen atom ) we resign of a relativistic description and combine 
quantum field theory with relativistic quantum mechanics. It seems that the 
conventional QFT is unable to treat any space-time description of relativistic 
particles. In particular, the limit h — ► of quantum field theory exists Q but 
it describes the classical field theory rather than the classical mechanics. Nev- 
ertheless, we would still be interested e.g. in the probability density of finding a 
relativistic electron close to the nucleus and a classical limit of this probability. 

In this paper we consider a non-linear relativistic quantum wave mechanics 
which is inseparably bound with the quantum field theory. The quantum effects 
are achieved by an addition of a noise to the relativistic field wave equation. It 
is shown that the expectation values of the random part of the field after an 
average over the proper time coincide with the time-ordered vacuum expectation 
values of the conventional QFT. 



II Relativistic wave equation 

We assume that the evolution of the wave function tj; T is determined by its 
initial value ip. In a classical limit we can imagine tp as a wave packet ^(a;) 
concentrated on a certain space-time point. We consider a relativistic particle in 
an electromagnetic field A. We suggest the following analogue of the Schrodinger 
equation (see earlier papers on such wave equations f 



*nd T i> = ^s^i-iM* + -^{-ihdv + X -A v )i> = -^a^ (i) 

where c is the velocity of light and <? M1/ = (1,1,1,-1) is the Minkowski metric ( 
we shall also use the notation R d for the Minkowski space suggesting an arbitrary 
dimension d of this space). If we write ( where W T may be complex ) 

4>r = exp(^W T ) (2) 



Then, from eq.(l) it follows 

drW T + + - c A,)(d»W T + \a») - ^<pW T + \d„A" ) = i) C-i) 

In a formal limit S-*0we obtain the Hamilton-Jacobi equation 



d r Wr + 2^(W- + \A^W T + ^4") = (4) 

Eq.(4) does not coincide with the conventional Hamilton-Jacobi equation in 
classical relativistic dynamics JtJ ^ which has no d T W term. Nevertheless, eq.(4) 
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can be derived from classical mechanics. Let us recall that if the relativistic 
Lagrangian is chosen in the form invariant under the reparametrization £(7) — > 
x (f(l)) then the canonical Hamiltonian 

H-^(p, + iW + ^) (5) 

is identically equal to zero. This constraint H=0 || generates correct equations 
of motion if the time parameter is interpreted as the proper time. If from 
the beginning we choose 7 as the proper time then H ^ 0. In such a case 
we may pose the problem of a canonical change of coordinates ( determined 
by the generating function W ) such that in the new coordinates H — > H + 
d T W — 0. The generating function W (hence also the classical dynamics) is 
defined by the solution of the Hamilton- Jacobi equation (4). Eq.(l) can be 
considered as a quantization of eq. (4) . In the standard quantization scheme of 
constrained systems || one argues that the quantum theory should be invariant 
under the choice of the parameter 7 ( reparametrization invariance). Hence, the 
dependence on this parameter is gauged away and what remains is the Klein- 
Gordon equation Hip = 0. However, in our interpretation the proper time has a 
physical meaning. Hence, we make this preferred choice of 7 in the Lagrangian. 
In such a case the canonical Hamiltonian (5) is different from zero. It again 
generates the correct equations of motion. Then, the conventional quantiztion 
scheme leads to eq.(l) rather than to the Klein-Gordon equation. We show 
below that conversely the classical dynamics (4) is determined as a limit Ti — ► 
of the quantum dynamics (1) if r is identified with the classical proper time. 

For further purposes we write the metric tensor in eq.(l) in terms of complex 
vierbeins v 

i«r = « (6) 

We can take as a solution of eq.(6) 
if k=l,2,3 and 



where 



The Feynman integral [Q supplies a simple intuitive way of proving the clas- 
sical limit. We assume that all functions we deal with are analytic. Then, we 
can express the solution of eq.(l) by the following rigorous form of the Feyn- 
man integral ( r > 0, see refs. Jl(]][[Ll] for a more precise formulation and all 
assumptions ) 

tp T (x) — E exp(-^- dsA fl (x + avb s )av^dbg)^(x + avb T ) (7) 



\~k fa dsA fi( x + <rvb s )avgdb ( *)ip(x + avb T ) 
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where b is the Brownian motion i.e. the real Gaussian process with values in 
R d and the covariance 

E[b a (s)b c (T)} =8 ac min(s,T) 

In cq.(7) 

a = 



h_ 

M 

We write eq.(l) and eq.(7) still in another (equivalent) form. Let W T be the 
solution of the Hamilton- Jacobi equation (4) with the initial condition W . We 
express the wave function ip as a product 

^ T = ex P CjW T )<S> T (8) 

If ip T is a solution of eq.(l) then <I> T fulfills the equation ( assuming the Lorentz 
gauge d^A^ = ) 

a T $ r = - l -M T = A n$T _ ±(d»w T + - c A^ T 1^ T $ T (9) 



where 

□ = fTd^ 

The equivalent form of the Feynman formula follows from eq.(9). So, if g M (r) is 
the solution of the equation ( < s < r ) 

% - —^( d n W ( T - s ' + \M9( a ))) da + ™»adb a (s) (10) 

where W(t) is the solution of the Hamilton- Jacobi equation (4), then the solu- 
tion of eq.(l) with the initial condition ip = exp(j-W)& reads 

W(x) =exp(jrW(T,x))E[exp(- J^aW(r - s,q(s))ds)^(q(T))\ (11) 
In the limit % — > of the stochastic process we obtain the flow (here < s < r) 

Till 0(/i) terms we have in the semi-classical approximation 

(13) 
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If we differentiate eq.(12) once more over s and make use of the Hamilton-Jacobi 
equation (4) then we conclude that £ fulfills the equation 



where 



with the boundary conditions 

£(s)u=o = x 



ds \s=t M c 



(d^W(x) + -A^x)) (15) 



We could also obtain eq.(14) directly from the Feynman formula (7). For this 
purpose we make a shift of variables b a — > b a + f a . Then, for any functional \ 
of the Brownian motion ( b(s), s < r ) we have ( the Cameron Martin formula 

H) 

E[ X {b)\ =E[exp( K - £ f a db a ^-^-ds) X (b + /)] (16) 

We apply the shift (16) to the Feynman formula (7) with ip = exp(j-W)Q and 
av f = £ — x , where £ is the solution of the equation (14) with the boundary 
conditions (15). Then, expanding in a we obtain the formula (13) till O[o). 

Summarizing the result (14) derived either from the Feynman formulas (7) 
and (11) or directly from the Schrodinger equation in the form (9) we can say 
that the limit Ti — > exists if and only if there exists £ such that the equations 
of motion (14) are satisfied. As an immediate but important consequence of the 
equations of motion (14) we obtain 

d% de _ 

ds 2 ds 

Hence, the square of the covariant velocity is time-independent. From the 
boundary condition (15) we obtain that 



f f = ^(Wx) + l Mx ) W W(x) + \ 



If in the Hamilton-Jacobi equation (4) 



W, - (17) 



Then 

ds ds 



c? ( 18) 
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Eq.(18) means that s ( hence also r in the Schrodinger equation) is the proper 
time. Moreover, this interpretation must be preserved in higher orders of h 
because the subsequent quasiclassical expansion is determined by the first or- 
der, ft is well-known from classical relativistic dynamics that the equations of 
motion (14) with a time parameter s have a correct form only if s coincides 
with the proper time. If the dependence of W T on r in eq.(8) is determined 
by eq.(17) then exp(j-W T (x)) also solves the Klein-Gordon equation ( in such a 
case we regain the conventional relativistic wave equation ) . On the other hand 
tpT ~ exp{ lT ^ ) means that the wave function oscillates rapidly with a very 
small period of oscillations . It is known from classical mechanics || that such 
additional rapid oscillations have little effect on the mean behavior in a slowly 
varying potential. We suggest that even if the condition (17) is not satisfied 
then the evolution in the proper time describes rapid oscillations which have 
little effect on the evolution in the Minkowski time for most physically relevant 
potentials A^. Hence, after an averaging over the proper time we obtain the 
conventional quantum theory. 

For an understanding of eq.(l) it is important to see its non-relativistic limit 
( the problem is discussed in a different way in ref. jl3| ) . The non-relativistic 
energy e is related to the relativistic energy p , mass M and the momentum p 
by the formula 

e = c(p 2 + M 2 c 2 )i - Mc 2 = + o(-) 



2M c 



Let 



tp T = exp( lMc — —)ip T 



Then, in the limit c->cowe obtain ( we write A — (A, V) i.e.Ao — V) 

ihd T ip T = -ihd t ip T + -r^i-itN + -A) Vt + V(x, t)ip r = Kip T (19) 
2M c 

Strictly speaking we should have omitted the -^A term in the limit c — ► oo 
but we keep it in order to be in agreement with the conventional procedure. If 
the potentials A and V are t-independent then we can express eq.(19) as the 
Schrodinger equation with a new time t = i + r. In such a case r is just a global 
shift of time in the non-relativistic quantum mechanics. 

If the potentials are time-dependent then we obtain Howland's description 
[ pT[ of the time evolution in time-dependent potentials . In such a case t is 
treated as a coordinate on an equal footing with x ( the eigenvalues of K are 
called quasienergies and constitute a standard tool in an investigation of time- 
dependent systems). 

The relation between the r-evolution and the conventional Schrodinger evo- 
lution U sc h is well-known [ful 

exp(—--K)ip(t, x) = U sch (t, t - T)i/>(t - r, x) (20) 
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In particular, it follows that the scattering theories in terms of exp(-^-K) and 
U sc h are equivalent. 

Ill Quantum free fields and free particles 

Let us consider a time evolution of the wave function with no interaction 

h 2 

lhd T <P T = -2^-°0r (21) 

An easy computation shows that the wave packet 

<j>{x) = j dpcj){p)exp{^p li x^) 

evolves into 

dp<t>{p)exp(-^-^)exp(-p^) (22) 
If 4>{p) is regular in h and its support is concentrated at p c e.g. 

4>(p) = expi^ - ^-(p - p c ) 2 - ^-(p +p 0c ) 2 ) 



then we can conclude that 



0,(3) « (23) 



Eq.(23) correctly describes the evolution of a relativistic wave packet if r is 
interpreted as the proper time. 

We can describe a wave packet of k-particles by a generalization of eq.(22) 

4>(x(l),....,x(k)) = J dp(l)....dp(fc)0(p(l),...,p(fc))eay(2^O>''(j')) 
Its time evolution is determined by a generalization of eq.(21) 

ihd T <Kx(l), ....,x(k)) = ~ £ z(fc)) (24) 

The solution of eq.(24) shows an independent free evolution of each particle 

Pc(j)r 



Xr(j) =x(j) + 



M 



By a free quantum field we understand an object which can describe any number 
of free particles as excitations of the field at any point of the space-time. By its 
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physical meaning the quantum field of a large intensity is supposed to behave 
as its classical counterpart ( keep in mind the example of the electromagnetic 
field ) . It should also describe fluctuations corresponding to the Heisenberg 
uncertainty principle. We represent these fluctuations by the Gaussian field B 
with the covariance 

E[B(t, x)B(s, y)] = mm(r, s)5(x - y) (25) 

We quantize the wave equation (21) by adding the noise ( when M — > oo then 
<j> is the Brownian motion) 

#(t, x) = ^ d 4>(t, x ) dr + V2hdB(r, x) (26) 

The solution of eq.(26) is a sum of two pieces ( see the subsequent section): the 
first part is the wave function (21) and the second one is a noise. We describe the 
quantum field by means of the random part relating the correlation functions 
of the random field to time-ordered products of non-commutative operators. 

The solution of eq. (26) depends on the proper time parameter r. We sug- 
gest that correlations which are observed in experiments result from an average 
over the rapid oscillations in the proper time. We shall show that time-ordered 
vacuum expectation values of the conventional QFT of massless particles coin- 
cide with an average over the proper time 

<4>{x l )....<t>{x k )>=lim T ^ 00 {T -to)- 1 ( dTE[<j>{T,x-i) (f>(T,x k )] (27) 

J T 

For massive particles with mass M we should add the term TpMcfr on the r.h.s. 
of eq.(26) i.e. -> ^{U - M 2 ) (in cq.(l) M fulfills the* role of the mass 

only in the non-relativistic limit). 

There is an alternative to the stochastic quantization (26) of the wave equa- 
tion (21). We could simply treat eq.(21) as an operator equation in the Fock 
space. Then , the conventional quantum free field is a particular T-independent 
solution of eq.(21). We could add an interaction and continue the proper time 
quantization in the operator formalism. We outline such an approach in the 
Appendix. 

IV Mathematical aspects of the linear stochas- 
tic equation 

We consider a general form of the linear stochastic differential equation 

d(f) T = -iA(f> T dT + V2hdB T (28) 
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Eq.(28) is treated as an equation in the Gelfand triple |H| 

S{R d ) C L 2 {R d ) C S'{R d ) 

where S'{R d ) is the Schwartz space of tempered distributions. A is a real self- 
adjoint operator in L 2 (R d ). Hence, there exists the unitary group U 

U(t) = exp(-irA) = cos(tA) - i sin(i\4) (29) 

The solution of eq. (28) with the initial condition at r reads 



V = U(t - T O )0 + 



V2h [ U(T-s)dB s (30) 

J Tn 



It is common to work with real processes. For this purpose let us decompose 
the complex field = 0i + i<p2 into its real and imaginary parts 



01 (t) = cos(A(t - tq))^ + sin(A{r - T a ))(j)2 + V2ti / cos(A{t - s))dB s (31) 
and 

02 (r) =cos(A(T-T ))(l>2-sin(A(T-To))<j>i-V2h sin(A(r - s))dB s (32) 



The solutions 0i(r) and 02 (r) determine the transition function ( from the 
initial point (0i, 02) to a set r C S' ) of the stochastic process 

P(r , 0i, 2 ; r, r) = P((0i(r), 2 (r)) S T) (33) 

The stochastic process T determines a solution of the differential equation 

<9 T $ T = h 2 Tr(D 2 <& T ) - i{A(j), D^> T ) (34) 

where 

(£>$(0), /) = Um^oe- 1 ($(0 + ef) - $(0)) 

is the Frechet derivative ; the second order derivative D 2 $ is an operator whose 
trace gives the Laplacian in an infinite number of dimensions ( for a diffusion 
in infinite dimensional spaces see Jl(|). In physicists' notation 

(A<j),D<f>) = [ dx{A<j))(x) " <J) 



and 

5 2 <& 



54>{x) 

Tr{D 2 <$>) = J dx- 



S(/)(x)S(/)(x) 
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The solution of the differential equation (34) with the initial condition $ reads 

$ T (fa = (35) 

where T (</») is the solution of eq.(28) with the real initial condition <f) at r . 
We can also express the solution by means of the transition probability (33) 

M<£) = J P^faQ^dfa^dfa^^+ifa^ (36) 

P(t, fa T) also solves eq.(34) with a <5-type initial condition . The transition 
probability P is defined by the Gaussian measure which is uniquely determined 
by its mean and covariance. The mean and the covariance can be calculated 
explicitly from eqs.(31)-(32). It can be seen from these equations that the 
average (27) over r does not exist separately for fa and fa. We show however 
that the average (27) does exist for the complex field <f> = fa + ifa. First of all, 
concerning the mean value 



T 

dr(f, exp{-iA{T - t q )) fa = i{{exp{-i{T - t )A) - l)^" 1 /, <t>) (37) 



Hence, the limit (27) of the mean value is equal to zero if A is invertible on f 
i.e. if / € Range(A). On the other hand if we consider A = — ^yD and the 
initial condition (f> satisfies the free wave equation 

□0 = (38) 

then such an initial wave function gives a non-trivial contribution to the average 
(27). We shall discuss this term later. Now, consider the covariance 

E[{{fa,f)-E[{fa,f)]){{fa,,f')-E[{fa,, /')])] = 

f (/, A- 1 (exp(-iA\r - r'\) - exp(-i(r + r' - 2r )A)) /') 

The average (27) is 

Umr^UT - ro)- 1 £ drE[((fa,f) - E[(fa, /)])(( fa, f) - E[{fa,f)])\ = 
^(f,A-\r) - ^limr^ooCT - ro)- 1 ^- 1 /, (exp(-2i(T - t q )A) - l)^" 1 /') 

(40) 

The limit of the second term on the r.h.s. is equal to zero under the assumption 
that A^ 1 f and A^ 1 f exist. So, we define A^ 1 first on a restricted set of 
functions ( with no support on the mass shell). However, it is not sufficient to 
define quantum fields only on test functions f with no support on the mass-shell ( 
p 2 = or more general p 2 = —M 2 for a massive particle when A = —h(0 — M 2 )). 
We have to specify the correlation functions also on the mass shell. In the 
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distribution theory Jig ] this means an extension of a linear functional to the 
whole of S. Defining the limit T — > oo is equivalent to a particular extension . 
There is a natural definition of this extension resulting from the solution (30) 
of eq.(28) and related to the treatment of indefinite integrals and especially 
oscillatory integrals. We make the replacement (e > ) 

A->A-ie (41) 

Then, the distribution 

□ _1 (x,y) = lim e ^oC2TTy d J d P exp(-ip(x-y)){-p 2 -it)- 1 = A F (x-y) (42) 

coincides with Feynman's causal function which is equal to the time-ordered 
vacuum expectation value of the real scalar free field 

<O\T(0(x)0(y))\O>=ihA F (x-y) (43) 

Note that if (x — y) 2 ^ in eq.(42) then Ap is a regular function. This property 
remains true for the kernel JCt{x — y) of the operator 

Kt = (exp(-2i(T- t )A) - l)yT 2 

( eq.(40) but now with A = — . Hence, if (x — y) 2 ^ then the limit 

T — > oo in eq.(40) holds true not only in a distributional sense but also for 
the operator kernels ( in particular K,t(x — y) tends to zero for every x and 
y if (x — y) 2 7^ 0). This stronger mode of convergence of the r-averages may 
be important for a convergence in a model with an interaction when there is a 
non-trivial renormalization and the distributional convergence is not sufficient. 
We can generalize this result to an average of an arbitrary number of fields 

< ((4> r , h) - EU T , /i)])((0 T , f 2 ) - E[(4, T , / 2 )]).... 

((0 T) f 2n ) - EUr, hn)]) >= h 2n E P air S U U ,k) Uh^- X fk) 

where the sum is over the product of all pairs (j, k) in agreement with the 
Gaussian integral combinatorics ( the expectation value (44) is equal to zero if 
the number of fields is odd). With A — —Ti(U — M 2 ) eq.(44) can be proved by 
means of explicit computations ( as in eqs.(39)-(40) ) through an application of 
the Fourier transform. 

The limit (27) exists on test functions such that A^ 1 f is well-defined i.e./ 6 
Range(A). We could instead of eq.(28) consider the equation for A(j> 

d{A(j>) = -iA{A(j>)dr + V2hAdB 

Then, we would not have the problem of an inverse in eq.(44) ( we just let 
fj — > Afj in eq.(44) ). We would have obtained the correlation functions of 
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A<fi = X- The ambiguity in the inversion problem arises if we wish to express <fi by 
X ( this problem can be considered as another linear stochastic equation) . Such 
an equation in S' poses the problem |l5] of an extension of a linear functional 
defined on a subspace of S to the whole of S. 

In this section we have discussed only the random field corresponding to the 
quantum real scalar field. If the quantum field has more components then we 
proportionally increase the number of components of the random field (as well 
as the number of Brownian motions in eq.(26)). So, for example we treat the 
complex scalar field as a real doublet {4>i,4>2) ■ The stochastic equation for a 
free electromagnetic field depends on whether we add the gauge fixing terms to 
the Lagrangian or not . Without any gauge fixing terms it reads 

d T Af, = hid^F^dr + V2hdB^ (45) 

If an external current is added to the Lagrangian then the stochastic equation 
takes the form 

d T A f _ l = ihd v F^dr - ihJ^dr + V^hdB^ (46) 

We would not obtain a finite limit T — -> oo of the correlation functions (27) 
of A without any gauge fixing. However, we may restrict ourselves to gauge 
invariant observables e.g. to . Then , we can easily rewrite eq.(46) in a 
gauge invariant form as an equation for F^ v 

d T F a ^ = ifid" (daFftu - d^F au )dT - iTi(d a J^ - d^J^dr + y/2h(d a dB M - d^dBa) 

We can still simplify the problem of solving this equation if we take a divergence 
of both sides 

d T d»F afM = -ihaF af ,dT-ihd' 1 (d a J^d f ,J a )dT+V2hd f "{d a dB f ,-d f ,dB a ) (47) 

Now, we can apply the formulas (29)-(30) where A = hD. The exponential 
of this operator and its kernel are well-known. We obtain a formula express- 
ing the T-averages of d v F av by a Gaussian field Xa(x) with known correlation 
functions. The gauge problem appears as the non-uniqueness of the potential 
A solving an equation d v F av = Xa ■ We discuss here these elementary aspects 
of linear stochastic equations because they will appear in a more complex form 
in quantum gravity discussed briefly at the end of this paper. 

Eq.(45) without the noise is treated as a wave equation for the photon. Not 
surprisingly the solution is complex and 

F^(t, x)F jk (r, ar) + *bfc(r, x)F ak (r, ar) (48) 

can be interpreted as the probability density (unnormalized) of finding a pho- 
ton in the space-time point x measured at the proper time r. Note that eq.(48) 
gives a generalization of the conventional statistical interpretation of the elec- 
tromagnetic field when we consider r-independent (real) solutions ( then eq.(48) 
defines the electromagnetic energy density ). 
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V General Lagrangians 

We consider now an interaction among relativistic fields <j) a described by a gen- 
eral Lagrangian C(4> a ). Its action integral is denoted L((f> a ) . The r-independent 
(no noise) solutions of the wave equations for <f) a should coincide with the clas- 
sical non-linear waves. A proper generalization of stochastic equations of sec. 4 
reads 

# (r, x) = ih—— -dr + V2hdB a (r, x) (49) 

6(f>a(T,X) 

where 

E[B a {r,x)B c (s,y)} = min(T, s)S ac 6(x - y) (50) 
We write eq.(49) in a symbolic form 

dej) = -iA(t>dr - igG(4>)dT + V2hdB (51) 

Let ( as in eq.(35) ) <& T (<A) = E[$>((j> T (</>))] then it follows from a general theory 
of stochastic equations [[16] that $ T (0) is a solution of the equation 

<9 r $ r = g<S> T = h 2 Tr(D 2 <5> T ) - i(A(f> + gG{<j>), £>$ T ) (52) 

We can solve eq.(52) by means of an expansion in g. 

oo 

$ r = 5"*r n) (53) 

n=0 

Then 

d T ^ = H 2 Tr(D 2 $W) - i(A(f>,D^) - i(G((b),D^-^) (54) 

The solution of this equation can be expressed by means of the transition func- 
tion (33) 

& T n \cb 1 +ic/ )2 ) = ~i [ P(r Q ,c/ )1 ,cl )2 ;T,d ( f>' 1 ,dcl > ' 2 )(G(cl > '),D^ T n ~ 1 '>) (55) 



where cf>' = (/)[ + i(j)' 2 . In this way we obtain a perturbative solution of eq.(52) ( 
the convergence of the series for interactions G(4>) of physical interest remains 
an open problem) . 

In the Appendix we outline a proper time quantization of interacting quan- 
tum fields in the Fock space. We show in the lowest order of the perturbation 
theory that such a quantization leads to the same results as the stochastic 
quantization (49). We supply also some non-perturbative arguments to this 
conjecture. However, we think that the stochastic approach although equiva- 
lent to the operator one avoids many difficulties related to non-commutativity. 
Moreover, it is useful for computations; in particular, for numerical simulations 
( see refs.[|l^ |l^] for numerical simulations of complex Langevin equations of 
the form (49) in the context of the stochastic quantization scheme of refs. [|l9| 

[H ) . 
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VI An average over the proper time 

When we average over the proper time r ( as in eq.(27)) then we obtain a linear 
functional / on a set of functions of fields 

/($) =< $ > (56) 

We are interested in obtaining equations which could determine the averaged 
values. For real processes and real diffusion equations the averaged value /($) 
can be expressed by a measure called an invariant measure. Let us recall this 
definition |Q. The solution <& r (0) = E[$(<j> T {(/)))] defines a semigroup & T (4>) = 
(P T $)(</>) ■ We say that v is an invariant measure if for a dense set of functions 

di/(0)(P T *)($ = / MflHti (57) 

So, we could say that computing the expectation values with respect to the 
noise B is equivalent to a functional integral with respect to v. 
Differentiating eq.(57) over time at r = we obtain 

J dv{4>)g${<j>) = J d{g*v)$ = (58) 

where the adjoint operator Q* is defined by the duality in the space of linear 
functionals. 

We shall show that the averaged value (27) is formally expressed by the 
Feynman integral. In fact, the averaging over the proper time could have been 
considered as a rigorous definition of the Feynman integral. We have already 
calculated the average values as the limit T — > oo of the correlations of the free 
fields (44). On the other hand a formal Feynman integral gives the same result 

d0esp(iL' ree fo))0(xi) <j>{x 2n ) =< <f>( Xl ) ct>{x 2n ) > /, ' ee 



here 



We shall show that in each order of the perturbation expansion the correlation 
functions defined by an average over the proper time fulfil the same differential 
equations as the ones resulting from a perturbative calculation of the Feynman 
integral ( the statement has been shown by means of a formal argument con- 
cerning distributional limits of oscillatory integrals in refs. (2lj|^|; the method 
applied here has been introduced first in the context of the stochastic quantiza- 
tion in ref.^3|). We neglect problems of the ultraviolet divergencies. In order to 
avoid the singularities we should regularize the noise B. Then, we prove that the 
r-averaged stochastic correlation functions coincide with the regularized corre- 
lation functions of the conventional QFT. In two-dimensional QFT we need no 



Stochastic quantum Held 



15 



ultraviolet regularization , hence the proof applies directly ( but because of the 
infrared problems we add the mass □ — > □ — M 2 ). 

As a first step in the proof we note that if the average over the proper time 
exists then from eq.(52) 

= UrriT^iT - to)- 1 f*dTd T E[*(<f> T (<f>))] = m 
= lim T ^ oo (T-T )- 1 J^Tg 4 ,E[$(<f> T (<l>))}=g <p <<S>(<i>T(<i>)) > 

If the average in the last line of eq.(59) is expressed by a (complex) measure v 
then eq.(57) defines an equation for this measure 

Q*v{<t>) = 

D(h 2 D + iA<j) + igG(<j)))v = (60) 
If A = -hD and 

G{4>) = hDV{(f>) (61) 
then the solution of cq.(60) as a formal complex Feynman measure reads 



or 



1 

dv((f) = d(j)exp(—L((f>)) 



where 



L{<j>) = J dxC{<j>(x)) = - J dxfynj + gV(4>{x)j) (62) 



VII Scalar quantum electrodynamics 

We discuss in more detail the scalar electrodynamics with {(jxj)) 2 interaction ( this 
interaction is needed for stability in less than four dimensions and additionally 
for renormalizability in four dimensions). We treat a complex scalar field as a 
doublet <p a (a — 1,2). Then eqs.(49) read ( in the Feynman gauge; e denotes 
the electric charge) 

d(j) a = ihO(j) a dT + igh{4>l + 4>l)<f) a dT - iet^A^d^dT - j-e 2 A^dr + V?hdB a 
= iha A( j) a dT + igh{(j}\ + 4>l)<j)adT + V2hdB a 

(63) 

dAy, = ihUA^dr + iheeabfibd^adr + V^hdB^ (64) 
We write eq.(64) in an integral form specifying the initial condition 

A^{t,x) = exp(ih(r - t )O)A^ + ieh J^ o exp(ih(r - s)U)e ab (f> b d ll (t>ads+ 
+V2hf ra ex P {ih{T - S )D)dB^( s ) = A<* + AQ 



Stochastic quantum field 



16 



where we denoted by Ajjj the noise-independent part of A^. Then, eq.(63) is 
rewritten in an integral form 

4>(t, x) — exp(ih(r — T )U A )(j) + igh J*J~ exp(ihd A (T — s)){4>\ + (/>|)(/>(s)ds+ 
+V2?i Stq exp(iha A ( T - sfjdB{s) 

(65) 

It is understood that <j> in this equation is a two-dimensional vector and A^ is 
a matrix with matrix elements (A^) a b — e a bA^ . If we insert as the initial 
conditions in eqs.(63)-(64) then as follows from sees. 6 and 8 after averaging over 
t we obtain the standard scalar QED |^4|. In particular, the formula (74) of the 
subsequent section can be applied for a computation of the scattering ampli- 
tudes. If the initial conditions in eq.(63) and A cl in eq.(64) are different from 
zero then our formalism goes beyond the conventional QFT. We restrict our- 
selves here to a discussion of the approximation with e = in eq. (64) whereas in 
eq.(63) we neglect the noise and the terms non-linear in <f>. Then, eq.(63) is just 
the relativistic wave equation (1) with the electromagnetic field A which is the 
sum 

A d + A Q In guch 

a case the solution (65) of eq.(63) is again expressed by 
eq.(7) and eq.(ll). We can compute the contribution of the quantum electro- 
magnetic fluctuations upon the meson scattering processes as well as the effect 
of these fluctuations on the energy levels of mesonic atoms ( the Lamb shift) . 
In this approximation to QFT we have a well-defined limit h — > of QFT when 
the first term on the r.h.s. of eq.(65) describes a localized trajectory of eq.(13). 
The method of calculations in the framework of the proper time stochastic field 
theory may be considered as a rigorous version of the argument of Welton p5f . 
We consider the approximation (10)-(13). The solution of the Hamilton-Jacobi 
equation (4) can be expressed by the solution of the Newton equation (14) 



W(r,x) = W(ar,x))+ (^^- + -A^)d» 



2 ds ds c 

We compute perturbatively the effect of A® upon the trajectory ( we performed 
similar computations in ref.pq]) . We write 



where 



and 



ds 2 W« ) ds 



M^ff- = (Frtf*+t™)-F*tf*))^+F^e*+ft^ (66) 

We can compute from eq.(66) perturbatively in powers of A® ( with 

A cl = (0,0,0,--^) we obtain Bethe's p7j and Welton's 5(x) correction to the 
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potential). Then, the evolution of the relativistic wave function in the electro- 
magnetic field A cl + A® is computed from the formula 

<f> T =« exp(ih(T — tq)O a )<P >> 

where << .. >> means the average over the electromagnetic field in the ap- 
proximation e = in eq.(64). 



VIII Large time behavior and scattering 

If the interaction is weak at large distances then we expect that the wave function 
4> T or the corresponding stochastic field behaves at large time as a solution of 
the equation with no interaction ( eq.(28) with A — > Ao, a differential operator 
with constant coefficients). Without noise such a behavior means 

(f) T w exp(-iTA )(f>in (67) 

when r — > — oo ( on a formal level such a behavior can be imposed if we rewrite 
eq.(49) in an integral form ) and 

(j) T k, exp(-iTA )(f> out 

when t — > +oo . 

In non-rclativistic quantum mechanics if the Schrodinger time evolution is 
determined by H then a local distortion of the behavior (scattering) is described 
by the operator 

S = Zimt-»oo,t -»-ooW(Mo) = /gg\ 
lim t - to0tto -,- 00 exp(^H t)exp(-^H(t - t ))exp(-j-H t ) 

If H — Hq = V then U satisfies the equation ( useful for perturbative calculations 

) 

ihd t U = V{t)U (69) 

where 

V(t) = exp(^H t)Vexp(-^H t) 

Similar formulas hold true for a relativistic quantum mechanics described by 
the Hamiltonian A 

Srel — HtTL t — >oo,To — > — oo^rel (^~; ^0 ) f7f\\ 

UmT^oz^^-oeexpiiAo^exp^iAiT - T ))exp(-iAoT ) 
If A — Ao = V then U re i satisfies the equation 

d T Urel = V{r)U rel (71) 
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where 

V(t) = exp(iAoT)Vexp(— iAqt) 

When c — > oo then ( as we have shown at the end of sec. 2 ) the evolution 
operator of eq.(l) tends to the evolution operator (19) of Howland. Then, it 
follows from the formula (20) ( see also |l4|]) that the scattering operator (70) 
of the evolution equation (20) in the limit c — > oo coincides with S of eq.(68). 
Hence, S re i — > S. 

We expect that in the scattering of particles of low energy the quantum 
held theory should give the same results as the quantum mechanics . We could 
establish this property in scalar QED of sec. 7 when the iterative solution A 
of eq.(64) starts from A cl such that UA cl = 0. Then, <f> T = exp(—ihTO Ac i)(f> 
plus the noise and higher order terms . If the noise and terms of higher order 
in the coupling constant are neglected then we obtain the time evolution and 
scattering of relativistic quantum mechanics. 

The quantum held theory is represented as a non-linear stochastic wave 
mechanics. The asymptotic free behavior (67) makes well sense in such a non- 
linear wave mechanics. The S-matrix can be defined as 

(f>out = S4> in (72) 

Its matrix elements are products of the in and out states. So, if the in- state 
f q (x) satisfies A f q = ( in cq.(51) we let A — >• A = -&(□ - M 2 )) then the 
matrix element is 

ifq,<t>out) = lim T _ 00 (f q , exp(-iTAo)<j> out ) « lim T ^oo / (f q ,d T 4> T )dT 

Jo 

We treat the multiparticle scattering as a scattering of independent particles 
connected only by the common noise. Then , a suggested generalization of the 
one particle matrix element to the multiparticle one could be 

<pi,...,Pn\S\qi,....,q m >=lim T ->oc7^ / drE\ JJ(/ Pi , 9 T (f) T )(f qk , d T <p, 

Jo 3 k 

. (73) 

If we replace d T 4> T by —iA^<f> T ( which follows by a formal differentiation of 
eq.(67)), average over r first ( before the differentiation contained in Aq = 
□ — M 2 ) then from the stochastic field theory we obtain the standard LSZ 
formula for the scattering matrix p4j 

<pi,...,p n \S\qi,....,q m >= J dxi dx n dyi....dy m f Pl (xi)....f Pn (x n )f qi {iji).. 

f qm {y m ){0 Xl — M 2 ) {U Vm -M 2 ) < <j} T {xi) <j) T (x n )(j) T (yi)..4 T {y m ) > 

(74) 
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IX Dissipative relativistic quantum mechanics 

Quantum theory of microscopic phenomena cannot be separated from a descrip- 
tion of the macroscopic world of measuring devices. An interaction of a quantum 
particle with a macroscopic environment is often described by a random wave 
function or (equivalently) by a density matrix 

p{x,y) = E[<f>{xj${y)] (75) 

If we require that the trace and positivity of p are preserved by the time evolution 
then we obtain the Lindblad equation p8| ( Lindblad equation for the relativistic 
quantum mechanics has been discussed earlier in ref.p9| ) 

d T p = -i[A, p\-\Y; R i R xp - \ E p R t R k + E R *p R t ( 76 ) 

fc k k 

where R describes a (phenomenological) dissipation. 

We can equivalently represent the dissipative dynamics by the Ito stochastic 
wave equation (i.e. a random perturbation of eq.(l)) which resembles eq.(26) 

d(f> = -iA4>dT + i Rk4>dB k - 2 E R k R k4>dr (77) 
fc fc 

where B k are independent complex Brownian motions 

E[B k B r ] = 

E\Ei(s,x)Bk(T,y)] = Skimin(T,s)d(x - y) 

If the operators R are Hermitian then eq.(77) can be expressed in a more com- 
pact form 

d<p = -iA<f)dT + i^2R k (f)o dB k (78) 

fc 

where the circle denotes the Stratonovitch differential |l2| and B k are indepen- 
dent real Brownian motions. When we solve eq.(77) and calculate the expecta- 
tion value (75) then we can see that the solution of eq.(77) determines a solution 
of the Lindblad equation (76). If the wave function <p has more components then 
these components add discrete indices to the density matrix. As an example 
, the random wave equation for the electromagnetic field in the gauge Aq = 
reads 

dA 3 = ihd F j0 dT - ihd k F jk dT + iR ljk A k dBi - \-Ri ko Ri km A m dr (79) 



We define 



p jk (x,y) = E[Aj(x)A k (y)} 
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Then the Lindblad equation takes the form 

drPjk ^[A : p\jk "^Rlrj Rlrm Pmk "^RlrmRlrk Pjrn ~\~ RlkmRljr Prm ($0) 

where 

Ajk;lm = Mj/(5fc m n — hdkmdjdl 

Eq. (80) can describe a time evolution of an ensemble of photons ( of various po- 
larizations corresponding to different indices of p) which undergoes a measure- 
ment R. In particular, a continuous observation of photon's position corresponds 
to Rijk = SjkXjdi ( where the vector a selects an orientation in space). 

Let us note that if we add a dissipation in the form (78) to the field equation 
(28) 

d(j> = -iA<j)dT + iRcj> odB + V2MB (81) 
and first calculate the expectation value over B 



p(x,y) = E[E B [<P(x)]E B [^y)} 

then we obtain again the density matrix p satisfying eq.(76) . However, we 
interprete eq.(81) in a different way. We consider eq.(81) as an equation for a 
quantum field interacting in a dissipative way with an environment. Then, both 
noises B and B arc treated on an equal footing i.e. we average over the noise 
at the end. We generalize eq.(81) to the non-linear quantum field theory (51) 
allowing R to be a real non-linear function of <j> 

d(f> = -iA<f>dr - igG{(f)dT + i^R{4>) o dB + V2hdB (82) 

Let (as in eq.(35) ) $ T (0) = E[®(<j> T (4>))] then it follows from cq.(82) that 
$ T (4>) is a solution of the equation 

2 

o> T $ T = G<S> T = fi 2 Tr{D 2 <S> T ) - ^Tr(R(<P)DR(<P)D<i> T ) - i(A<p + gG{(j>), D$ T ) 

(83) 

We average again over r and look for a measure v such that an average over 
v is equivalent to the r-average. We obtain a condition which is an analog of 
eq.(58) 

J dv{<(>)g$(cf>) = J d{<3*v)§ = (84) 
We rewrite eq.(84) as a differential equation for v 

2 2 

D (h 2 D - l-R^R^D + l-R{cj>)DR{<p) + iA<j> + igG{<j>j^ v = (85) 
Let as in eq.(60) A = —ho and 



G{<p) = hDV(<p) 



(86) 
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We express the solution of eq.(85) as a formal complex Feynman measure 



dv(<f)) = d(j)exp( — L^(4>)) 



where 



L 7 (0) = J dxCWx)) = J ' dip( - ^U<j> - gV(4>(x)) + \H<£)) (87) 

Then, eq.(85) for v can be rewritten as a linear equation for T 

fi 2 DT = i?(0) J R(0)DL 7 + ihR(4>)DR(4>) (88) 

Its solution is a complex function of <fi ■ The imaginary part of T describes a 
dissipation. In the lowest order in 7 we obtain 

The Feynman integral representing the averaged correlation functions (27) is an 
oscillatory integral describing the interference phenomena in quantum mechan- 
ics. The addition of a dissipation damps the interference leading to decoherence 
| j30| and a smooth equilibrium limit |jl| represented by the complex Gibbs mea- 
sure (87). 



X Quantum fields on a general curved manifold 

The proper time formalism is especially useful if fields are to be defined on a 
pseudoriemannian manifold M. . In such a case there is no candidate for a 
time as an evolution parameter ( the coordinate xq is only locally defined ). 
Moreover, the classical proper time (intrinsic time) is a function of the metric. 
Hence, in quantum gravity , when the metric becomes a dynamical variable, 
then the proper time acquires the same dynamical content as the coordinate 
x in non-relativistic quantum mechanics. We can easily generalize eq.(l) to a 
general pseudoriemannian manifold M ( we set c = 1 in this section) 

ihd T Tjj = ^9^(~^ + A^ + hg^T^i-ihdv + A v )^ (89) 

where g^ v is the Riemannian metric and T is the Christoffcl symbol of the 
Levi-Civita connection . Assume that W T is a solution of the Hamilton-Jacobi 
equation 

d T W T + J-g^^Wr + A^){d v W T + A,) = (90) 
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with the initial condition W. Let us consider eq.(89) with the initial condition 
</> = exp(j-W)& . Then, (f) T is a solution of eq.(89) if and only if $ T is the 
solution of the equation ( the Lorentz gauge for is assumed ) 

ihd T $> = -^n fl $ - Qg^PuWr + A^d^r - ^n ff $ T (91) 

where D g is the wave operator on the pseudoriemannian manifold M. 

In the formal limit ft — > of eq.(91) we obtain (j> T (x) ~ $(£(t)) where £ is the 
solution of the equation (0 < s < r) 

^ = -^m)(duW(r - s,t(s)) + A v (£(s))) (92) 

Differentiating eq.(92) once more and using the Hamilton- Jacobi equation (90) 
we obtain the equation 

^ + r ^^ = ^™ s) )^ (93) 

ds 2 + up ds ds M { ^ { )} ds ( ' 

Eq.(93) shows that the correct classical limit ( as a geodesic equation ) of a 
motion of the quantum particle on a general pseudoriemannian manifold results 
if and only if t is the proper time. The interpretation of r as the classical proper 
time remains true not only in the leading order in h but also in all subsequent 
terms because the leading order determines the subsequent interpretation of r. 

We can express the quantum corrections in terms of the Brownian motion. 
Let us define the complex vierbeins v^a 

V^aVva = i^v (94) 

Then, we consider a complex Markov process q as a solution of the set of 
covariant equations ( for the case of the Ricmannian manifold and imaginary 
time see Jl^| or |57j ; the real time is discussed in , it can be considered as a 
complexification of the diffusion ) 

d<f s = -^9^(qs)(d u W(r -s,q s ) + A v (q s ))ds + av%(q s ) o db a (95) 

d< + r^ p <°^ = o (96) 

Then, the solution of eq.(91) with the initial condition -0 = exp(j-W)^ reads 

i> r (x) = exp{^W T (x))E[exp(- — J a g W[r- «,g(«))ds)*(g(r))] (97) 
Eq.(97) can be applied for a semiclassical expansion in powers of h. 
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We can continue now (as in the earlier sections) with the quantization of the 
wave equation (89). So, for the free field the stochastic equation coincides with 
eq.(28) where 

A = -7iO g (98) 

In the computation of an average over the proper time (27) we need to define 
Dg 1 - At finite T we obtain a linear functional on a subset of test-functions 
S(A4). Taking the limit T — > oo is equivalent to a definition of an extension 
of this linear functional. The definition of the inverse is not unique but it is 
restricted by the requirements that the quantum field theory resulting from D" 1 
is causal and that the quantum fields should be defined in a Hilbert space ( with 
a non-negative scalar product). It remains unclear whether such a definition 
of 1 exists on an arbitrary manifold. This can be achieved on a globally 
hyperbolic manifold J|^] |53| |35| . In such a case there exists a complete set 
of solutions of the wave equation 

a gUj = (99) 

Then, we can define 

& { + ) {x,y) = Y J U 3 {x)u J (y) (100) 

3 

Eq.(lOO) determines a non-negative bilinear form. We can define by means 
of the standard Fock space methods the quantum annihilation <p^{x) ( as 
</>W = cijUj) and creation <f>^\x) parts of the field operator <fi — <fi( + > + </>("). 
Now, on a globally hyperbolic manifold there exists a choice of the xo-coordinate 
[|4| such that 

U-\x,y) = A F (x,y) = 9(x - y Q )A g +\x,y) + % - x )A^(y,x) (101) 

is independent of the choice of coordinates. Eq.(lOl) defines the Feynman prop- 
agator on a globally hyperbolic manifold . The definition (101) coincides with 
the one resulting from the e-prescription p6| (e > 0) 

/•OO 

□ ~ 1 = lim e ^o / idsexp(—isO g — es) (102) 
Jo 

The interacting field is defined by the stochastic equation (49) with A — —TiU g . 
After the propagator is defined there is no further difficulty (up to the ultraviolet 
problems) in defining perturbatively an interacting field by means of an iterative 
solution of the stochastic equation (49). 

A really ambitious program concerns a quantization of the metric g. The 
proper time approach supplies a useful method at the problem where other quan- 
tization methods encounter insurmountable difficulties. The Feynman integral 
is a formal tool . It can be defined either in the imaginary time or for complex 
coordinates |fl0|| . However, these methods fail for Einstein gravity which has the 
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action unbounded from below. The method of Euclidean stochastic quantiza- 
tion |59]]|4(| works well . However, it remains unclear whether the model can be 
continued back to the real time ( pseudoriemannian manifold) . An application 
of the stochastic quantization JlS| ] with a complex action and real time to the 
Einstein gravity has been suggested by Rumpf p8j . We interprete the fictitious 
time of Parisi and Wu as the proper time with a physical meaning. Such an 
interpretation should have experimental consequences for a time evolution of 
the graviton wave function ( in ref. [p~t| an experiment showing the interference 
in time is suggested). 

We need now the Brownian motion B which is a symmetric matrix whose 
entries constitute independent Brownian motions 

E[B ac (s, x)B fd (r, y)] = (6 af 5 cd + S ad S cf )min(T, s)5(x - y) 

Then, we introduce the infinite dimensional pseudoriemannian manifold Q{M.) 
of metric tensors g on M. with DeWitt's supermetric |[l2] on this manifold 

G a ^ = (det(-s))-* G^y* + ff "V a ) (103) 
We need a basis of complex frames on Q(M) (tetrads) fulfilling the relation 

If L is the action for the metric and the matter fields then the stochastic 
equation for g according to the prescriptions (49) and (95)-(96) reads ( see refs. 



7) ]39 ]| 40 for the case of infinite dimensional Riemannian manifolds ) 



c t 

dg^ = ih- dT + Ti£^ ac o dB ac (104) 

og^ v 

d £^.,ac + r[^ ) )(7p) £ 7P;aC O dg^ = (105) 

Eqs.(104)-(105) need a regularization of the Brownian motion B for a proper 
interpretation of the Christoffel symbol appearing in eq.(105) ( see (39|). Note 
that if L is a sum of the Einstein action and the action for the matter fields 
then eq.(104) reads 

dg» v = in(det(-g))i(R" u - ^g^R)dr - ihn{det{-g))^T> lv dr + h£^ ac o dB ac 

(106) 

where k is proportional to the gravitational constant and T^ u is the energy- 
momentum tensor for the matter fields. If we start from a linearized form of 
cq.(106) then the first difficulty which we encounter involves the zero modes of 
the linear part which cause a trouble with the large time limit. We can deal 
with the problem of zero modes cither by means of the stochastic gauge fixing 
[E3| or restricting ourselves to gauge invariant observables. We would like to 
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outline here the method introduced for the electromagnetic field at eq. (47) ( the 
Einstein tensor G^ v = (det(— g))* {R^ u - \g^ v R) is an analog of , it has a 
gauge invariance resembling the one of non-Abelian gauge fields p3|). So, for 
pure gravity ( T^ v = 0) we have 

d((det(-g))HR^ - \gTR)) = A ° dg a0 (107) 

2 ogapdg^v 



or more explicitly 



d({det{-g))^{R^ v - \g^R)^j = 

lh T^kz ((det(-g))i(R a P - \g a0 R))dr - i£ a ^ ac o dB a 



(108) 



Let g^iv — g\j,} + hfj,u where g^J is the constant metric tensor of the Minkowski 
space-time. Denote by G(h) the part of the Einstein tensor linear in h. Then 
eq.(108) reads 

dG» v (h) = 2ihUG^{h)dT + p{G{h)Y v dT + C» v ' ac {G(h)) o dB ac (109) 

where p and C are non- linear in G(h). In eq.(109) we have inverted the linear 
relation h — > G(h) i.e. we have chosen a particular solution /i M „ of the equation 
G^ u {h) — G^ u . We can solve eq.(109) iteratively ( we assume a coordinate 
invariant regularization of the ultraviolet divergencies J39| ) . Then , the proper 
time average (27) of G^ v (h) exists in a perturbation expansion with defined 
by /\p(x,y). We compute in this way the time-averaged correlation functions 
of G (h). These correlation functions determine a random variable x^v- Then, 
we can obtain the metric tensor as a particular solution of the linear equation 
G>„(/i) = Xfj,i> ■ Subsequently, we can express g^ v in terms of \nv ■ We expect 
that in spite of the arbitrariness (gauge dependence) of the relation between 
G(h) and h the coordinate independent variables e.g. R liUpl (g)R^ l ' pl (g) ( where 

g^v — + hsjui and R is the Riemannian curvature ) will not depend on the 
choice of the gauge ( as in the case of the stochastic gauge fixing). In general, 
we could start from a background metric gj±J instead of the Minkowski one. In 
such a case in addition to the processes of graviton creation and annihilation ( 
and mutual scattering) we could describe the time evolution and localization of 
a single graviton in interaction with other fields. 



XI Summary 

We have discussed a method of quantization which can be considered as a version 
of the Feynman integral. The method is related to the stochastic quantization 
of Parisi and Wu JljJ . We considered the Bose fields only. If we are interested 
in correlation functions of Fermi fields then we need stochastic equations with 
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Grassmann variables. Our main emphasis was on the proper time interpretation 
of the "fictitious time parameter". The proper time interpretation supplies a 
relativistically covariant quantization. It can give a meaning to the probability 
of detecting a particle inside a bounded region (sec. 7) . It is rather difficult to 
formulate such questions in the framework of the conventional QFT. Moreover, 
a proper time interpretation can have experimental consequences for time mea- 
surements |ffl|| . For quantum field theory on a manifold and quantum gravity it 
may be valuable to treat the background metric and gravitons at the same foot- 
ing. The conventional Euclidean quantization of Einstein gravity is impossible 
because the action is unbounded from below. The stochastic Euclidean quan- 
tization of gravity is feasible |^9||^]. However, in such an approach it remains 
unclear how to continue to the real physical time. The direct quantization in 
the real time avoids these problems. We have clarified some mathematical as- 
pects of the proper time stochastic quantization scheme. We have shown the 
convergence of the averaged values. The convergence rate of the average over 
the time interval [0,T] is usually proportional to T _1 . With purely Hamilto- 
nian dynamics it cannot be as fast as in the Euclidean version of the stochastic 
quantization where the convergence can be exponential. Nevertheless, in view 
of the difficulties of the conventional oscillatory Feynman integral the stochastic 
version can be a useful tool for numerical calculations ( stochastic simulations 
of complex stochastic equations are discussed in [{L7]|]l8|][fLl||). In the dissipa- 
tive quantum field theory of sec. 9 the convergence of the averaged correlation 
functions can again be exponential as in the Euclidean framework. 

XII Appendix 

We outline in this section an operator quantization in the proper time. The 
equation of motion reads 

SL 

d T <t>{T 1 x)=ih— , (Al) 

54>(t,x) 

We are looking for a solution of this equation with an initial condition at — oo 

lim T ^^ 00 (l) T (x) = 4> in (x) (A.2) 

where (□ — M 2 )4>i n =Q is the quantum scalar free field of mass M. We expect 
that as r — > oo 

ct> T (x) -> (/) ut(x) (A3) 
If this asymptotic behavior holds true then we can define 

4>out = S~ 1 cj) m S (A A) 
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We can develop a perturbative theory now. Assume the Lagrangian (62), 
then we can rewrite eq.(A.l) as an integral equation 

(j) T (x) = (j> in (x) - igh [ exp( - ih(U - M 2 )(r - s)\ V(<p s )ds (A.5) 

J —oo 

We are interested in computing the time-ordered products of quantum fields. 
We first solve eq.(A.5) perturbativcly ( fixing the initial condition at r — ► — oo) 

4>t = (pin - igh fr e«p( - iK u - M 2 )(t - s)jV(<j) m )ds 

-(igh) 2 X o expi^ - ih{r - s)(D - M 2 )^ V(4> in )ds (A.6) 

/; o exp( - ih(n - M 2 ){s - s'j)V{<j> in )ds' + .... 

We assume that V is defined by the Wick normal product. Then , the Wick 
theorem allows us to calculate explicitly < O\T((f> T (x)^) T '(x'))\0 > where the 
T-ordering is understood in xo rather than in r. Subsequently the average 
over the proper time can be calculated. We obtain a closed formula in terms of 
the Green's functions of exp{— is(U — M 2 )) and (□ — M 2 ) -1 . We must define 
□ [x, y) as Ap(x, y) if the time-ordered correlation functions are to agree with 
the conventional ones. We have checked the agreement with the conventional 
QFT only at the lowest order of g. We have obtained the conventional QFT 
because of the asymptotic condition (A. 2). Eq.(A.l) at g — has more solutions 
which in general are r-dependent. Starting with such solutions we would go 
beyond the conventional QFT. 

Another ( implicit) argument concerning the relation to the conventional 
QFT at each order of g comes from the identity 

= limR^ocjz J^dr < T(d T (4> T (xi) (j) T (x„))j >= UmR^^^ f^dr 

ELi <t( k M*i)--{(v - M 2 )^k) - gV'(Mxk))--4Ax n )) > 

(A7) 

We can use this equation in order to prove that if at g = the time-ordered 
correlation functions coincide with the conventional free ones then at each order 
of g the proper time average coincides with the one resulting from the solution 
of the operator equation 

(□- M 2 )0(x) = gV'Mx)) 
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